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Abstract 

The letter presents new field-theoretical approach to 2D passive scalar prob- 
lem. The Gaussian form of the distribution for the Lyapunov exponent is 

derived and its parameters are found explicitly. 
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DEFINITION OF THE MODEL 



The convection of a passive scalar (PS) in 2D by an external long-range velocity field is 
well-defined linear problem . To solve the problem is to find all correlation functions of 
the passive scalar field w(r; t) governed by the following equation 

Co + u a V a u = 0, (1) 

with u(r;t) being an external Eulerian long-range (q < L^ 1 ) velocity field and being a 
random external source localized in k space at fco = L~ x . The measure of averaging over 
the random source in the simplest case of white-noise source statistics can be chosen in the 
following form 

£™ <t>l =l / L {t)dt _ 2 r 

£>0exp( ) , (j) k = (2tt) J 0(r; t) exp()kr)dr, 

< Mt)M*) >s= 2k_* (* - t') = J\-8(t - t')5(k - k')S{k - k ), (2) 

without any loss of generality as it was shown in 0]. Here P 2 stands for the flow of the 
squared PS on the source. 

In order to eliminate a homogeneous sweeping it is convenient to pass to the locally 
comoving frame [[| expressing the Eq. (|I]) in terms of the Quasi-Lagrangian (QL) velocities 
related to the initial Eulerian ones by means of the formula 

u(r;t) = v(r - y < v(0;r)dr;t). (3) 

The eq.(|l|) takes correspondingly the form 

Co + K - <)V a u; = 0, , v = v(0; r). (4) 

Following the ideas of recent paper one can expand v a (r) — v a (0) = for the points 

satisfying r < L. Here a is the random in time traceless symmetric matrix. For the 
white-noise velocity field statistics (which is the only case considered in the present paper) 
the a elements measure has the form providing space-isotropy of the simultaneous pair 
correlator of the QL velocities 
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<a b ^ 



Va(t) = V 
To solve the resulting equation 



i b —a j 



VaVbexp(-S ),S = ± J(a 2 + b 2 )dt, (5) 



6j + a af} r p V a u = (f>, (6) 



one performs the substitution u(r;t) = ip(R;t), R a = W a/3 (t)r^ that allows to write down 
the solution of the eq. (|J) in the following form 

u(f;t) = ip(R(t);t) = f dr<p(W-\T, -oo)R(t);r) = f dr<f)(W{t,T)r;T), (7) 

J~ oo J— oo 

where the matrix W(t; r) can be sought for as satisfying the relation 

W t + Wa = 0, (8) 

with initial condition W(t,t) = 1. Due to Eq.(^) the all order correlation functions of the 
PS field may be rewritten in terms of the correlation functions of the source field 0. For 
example, the pair correlation function is 

< w(n; ti)w(r 2 ; t 2 ) >= I * dr x I dr 2 E(W(t 1 , n)ri - W(t 2 , r 2 )r 2 ; r x - r 2 ) a . (9) 

J —oo J —oo 

The explicit form of the simultaneous pair correlator ( which has to be isotropic) can be 
extracted from Eqs.© and (@) as 

< W (r i; t)u(v 2 ; t) >= P 2 J o °° dtJ ( lW{tM p~ r2)l ), (10) 

with Jq(x) being the Bessel function of zeroth order. 

So, following the calculations of the paper ||, we reduced our initial problem to the 
evaluation of some functions of the matrix W satisfying the Eq.(||). The statistics of matrix 
elements of W(t, 0) is defined by the ensemble of random real traceless symmetric matrices 
a, see Eq.(^). The authors of the paper have calculated <| W(t, 0)r |> by means of 
the direct expansion of the anti-chronological T— exponent that is the formal solution of 
the Eq.(||). It gave them possibility to estimate the simultaneous pair correlator (0) at 



small distances r\ — r2 as ^ ln(-). In addition, analyzing the set of many-point correlation 
functions of PS field Falkovich and Lebedev were able to demonstrate that the instantaneous 
statistics of such a field becomes more close to the gaussian one when passing down-scales. 

In the present paper we suggest a new nonperturbative field-theoretical approach to 
the calculation of the PS-field correlations. We show that the asymptotic behavior of the 
pair correlator Eq. ( [XQ|) is governed by the self-averaging (or "deterministic") quantity: 
the Lyapunov exponent. By means of our method we rederive, extend and exactly prove all 
abovementioned results by Falkovich and Lebedev. In particular, we find explicitly the value 
of the variance of the gaussian distribution characterizing the short range PS-field statistics. 



FUNCTIONAL REPRESENTATION FOR AVERAGED FUNCTION ALS OF W 



The matrix W(t) can be extracted from Eq.(|]) only in a form of the anti-chronological 
time-ordered exponent 

W{t,0) =f exp(- f dt'a(t')) , M0,0) = 1. (11) 
J o 

rather than in terms of some regular function of a. A similar problem - transformation of 
a time-ordered exponent of some linear combination of spin SU(2) operators (arising when 
one tries to write down an exact functional representation for the partition function for 
Quantum Heisenberg ferromagnet) has been solved by Kolokolov 0. The main idea of the 
method is to introduce such a new set of integration variables in the functional integral that 
Texp becomes some regular function when expressed in their terms. In the present context 
we use a new modification of the Ansatz proposed in |7[] expanding the a matrix over the 
spin 2x2 matrices as 



a = aa z + ba x . 



(12) 



First, we introduce a new basis of the spin algebra: 



f 



a„ 



a + =a z + io x 



V / 



a, - ia x 



-i -1 



;i3) 



that corresponds to the rotation of the quantization axis from the usual orientation (parallel 
to the z-axis) to the new one parallel to the y-axis. Instead of the fields a(t), b(t) we choose 
the new ones uo^ = (a ± ib)/2 transforming a and the integration measure Eq.(|5|) to the 
more compact form: 



2 r+oo 

a = <p~cr+ + y? + <x_ ; Va(t) = VyT exp(-S'o) , S = — / (f + ip~dt, (14) 

JJ Jo 

Let us now introduce the operator given in the explicit form as 



A(t,0) = exp f-o-_-0 + (O) 



exp 



-a a 



r(t')e 2 ti' p(t " )dt "dt' 



exp 



a y J p(t')dt' 



exp 



<7_<0 + (t) 



(15) 



where V ;± )P are some new dynamical fields. Using the commutations relations for the spin 
operators a±, o y it is easy to check that the operator A obeys the differential equation 

A t = A[-a+ip- + a-(A4)-(<4) + ) 2 + 2p^+ + + d y {-^~^ + p)}, (16) 

Using the fact that the first exponential factor in Eq. (p^5|) makes the operator A satisfying 
the condition A(0, 0) = 1 and comparing the last equation flTBD with the equivalent one (§) 
we find out the change of variables 



up =ijj , 

ip + = -ip+ + 4^"(^+) 2 ; 
p = 4ap~ip + , 



(17) 
(18) 



in case of which dynamical matrices W and A coincide. It allows us to obtain the explicit 
functional integral representation for any average written in terms of W by changing variables 
from ud ± to ip ± . The Jacobian of the map Eq.([H]) is 



(19) 



and it is essentially dependent on the kind of regularization and a type of conditions imposed 
on the field if) + . The matter is that the transformation Eq. (|i"7D contains the derivative of the 



field ip + with respect to time in its right-hand side. Therefore, it should be supplied with 
some initial or boundary conditions. In the papers |7| it has been shown that only imposing 
conditions of the initial kind it is possible to ensure the invertability of the map Eq. (|T7j) . 
Another point is that in the course of calculation it is necessary to average some functions 
of the operator A(T, 0) at fixed time moment T over the measure Eq. flT4]) . For a given T it 
is convenient to fix the final value of the field ip + 



^+(T) 



1 

2' 



(A 



\ / 



Then the operator A(T, 0) acting on the initial vector 
expressions: 

/ \ / 

1 -2^+(0) 
v i(l+2V+(0) j 



(20) 



produces the following simple 



R(T) = A(T, 0) 
R 2 (T) = -2^+(0)e 8 /o ^ 



(A I 

2 



\ / 



dt 



Here we exploited the isotropy condition and substituted R(0) 



(A 



vV 



(21) 
(22) 

for the initial value 



of the vector R(T) without any loss of generality. 

The regularization of the map Eq. (|l7|) (again as in the papers is determined by 

the physical argumentation. Here it stems from the translation in time invariance of the 
white- noise measure Eq.(|l4]). Indeed, the equality 



< (p + (T) [ T (p-(t)dt >=< <p~(T) [ T (p + (t)dt>, 



(23) 



leads to the extension of the definition of the step function 9{x) at x = such that 6(0) = 1/2. 



Thus, the discrete version of the map ([HI) = <^ ± (t n );n = l,...,M;e 
en; M — ► oo), compatible with the symmetry condition Eq . (|23|) is 



M 



0]t n 



(24) 



that gives the following expression for the Jacobian 



rT 

J = const exp(4 / ip + ip-dt'). (25) 
jo 

When calculating the Jacobian Eq . (p5|) the fields ip and ip^ were considered to be 
independent complex variables or, in other words, as different coordinate systems in the 
whole space C 2M of the fields' configurations. The conditions 

V + = (yO*, (26) 

being externally imposed on the model specify the surface £ in C 2M along which the differen- 
tial forms T>(p + f\ T>(p~ or f\ T>ip~ are integrated. For the coordinate sets (if) ) the equa- 
tion (|Tj]) for £ can be considered to be an implicit one. According to the Cauchy-Poincare 
theorem the integration surface can be deformed in an arbitrary way in the convergence 
domain provided we integrate an analytical function. There exists a continuous family of 
surfaces (homotopy) which includes both the surfaces £ and the "standard" one £': 

£' = {-^ + = (^r}. (27) 

The explicit expression for such a homotopy in a more general case can be found in the 
paper ||. Thus, we can replace the surface of integration £ by the standard one 

Substituting Eq.(|17|) into the measure Eq.(|) and using the expressions Eqs. ( |19p5"l) we 
obtain the following modification of the measure for averaging over ifj^ [[10] : 

r+oo r 2 • 



NV^expi-S^}) , St= [ 

Jo 



or , correspondingly, in the discrete form 



dt, (28) 



Si = £[-£(V£u - + "iUlfcDW + ^n + i) 2 " 2e0-(V£ + Vtf+i)]- (29) 

n 

It means that we reformulated the initial problem of the multiplicative random matrix pro- 
cess with the measure Eq. (j5j) to the multiplicative random scalar process with the measure 
given by Eq. (f2§|) . 
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GAUSSINESS OF PASSIVE SCALAR CORRELATIONS 



Another peculiarity of an arbitrary multiplicative random process is the Gaussian-like 
fluctuations of the Lyapunov exponent A(T) with an amplitude that decays like const /VT 
when T goes to infinity. Let us note that this fact has been proven both for the multiplicative 
random scalar process (see, e.g. JHJ) and for matrix one (see [fL2] .) The usual way to calculate 
the constant in front of 1/T 1 / 2 is to restore the distribution function of A(i) at arbitrary 
time from the set of moments (R 2n ) and then to extract the Gaussian-like peak around the 
mean Lyapunov exponent at large enough time (see for example the calculation of such a 
kind in ID [|L| and quasi ID [14| localization problem). We follow the same way in the 



present context as well. 

The Lyapunov exponent in our case is convenient to define as 

A,». (30) 
2T v ; 

In order to find the probability distribution function (PDF) of A let us firstly calculate all 
the even moments of R 

« #• >= siz^m^r^i, (31) 

< 1 >i 

where < ... >i stands for the average with respect to the measure Eq. (fffi). Passing to the 
new integration variables by means of the gauge transformations 

^ = X ± exp[ T 4 jf X + X~dt' ± (2m + l)D(T - t)], (32) 

that can be written in the discretized form as 

M 

^ = xt exp[-2 <X++i + Xf)xj + (2rn + l)eD{M - n + 1)] 

j=n 
M 

= Xn eX P[2 e (xt+l + Xj)Xj + e (Xn + Xt+l)Xn ~ ( 2m + l)eD(M - 7l) - 

j=n+l 

\e{2m + 1)D], (33) 



S 



we reduce the averaging in both numerator and denominator (the latter case corresponds to 
m — 0) of Eq.(|3T|) to the Gaussian- type ones 

< (-2^ + (0)) m e 8m /o T ^" d * > 1= (-2) m < (x + (0)) m > 2 e [{2m+2)m+1 z ]DT , < 1 > 1= e^, 

(34) 

Here (...)2 stands for the averaging with respect to the Gaussian measure 

9 r+oo 

N'Vx^expi-S^}) , S 2 = - / x + X~dt, (35) 

D Jo 

and one takes into account that the Jacobian of the map Eq.(|32"D is 

J[ X ] = exp[-2 [ T X + X~ + (m + \)D{T - t)]. (36) 
jo z 

The average < (x + (0)) m > 2 is equal to (— |) m due to condition x + (T) = — |. This result 
is easy to get shifting x + — > — \ + X + an d noticing that all the averages of x + are equal to 
zero. Thus, we arrive at the result 

R m (T) = e DT (2 m +2)m_ (37) 
Knowing R m one can extract Fourier-representation for PDF of R 2 in the following way: 

~ ~ {lk) m „ r+oo e ~DT/2 2 

and the PDF for R 2 is obtained after the Fourier inversion in the form 

—e~ ikz V k = / dx - exp[-— + -}5(R 2 - e x ). (39) 

Correspondingly, the PDF the Lyapunov exponent defined in Eq.(30) is explicitly given by 
the expression 



V(X) = ^e-^, (40) 

that means exact Gaussian fluctuations around the average value < A >= D with the 
variance ^ vanishing when T tends to infinity. In particular, A tends asymptotically to a 
non-random (" deterministic") quantity in agreement with general theory. 
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In order to pass from the statistics of A to the statistics of the PS-field we introduce the 
fluctuating quantity 



Q{j) = P2 J Q ~ dtJ o(je xt ), (41) 
whose distribution can be restored from the PDF of A. Indeed, at | C 1 we can cut the 



integral in the rhs of the Eq. at t = ln( -'' A //L ^ and estimate Q with logarithmic accuracy 



as 



«l)~i^§M)Mf). (42) 



It gives us the possibility to extract PDF of Q from the expression fl40|) . We obtain 



p 2 ,f^Hf) em 2P 2 g 



HQ) - , , T *M- D}. (43) 



Such a PDF has at r —>■ asymptotically the form of the Gaussian distribution for the 
quantity approaching the 5-functional form when j — > 00. 
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